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Abstract 

The scalar decay into gravitinos is studied in the presence of D-term SUSY breaking. 
The gravitinos produced by the decay of coherent oscillations of scalar fields, such as 
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moduli and inflatons, cause cosmological problems. We show the general formula for 
the partial decay rate of the scalar into two gravitinos. Applying it to a concrete model 
of D-term SUSY breaking, we find that D-term SUSY breaking can suppress gravitino 
production. 
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1 Introduction 



In models of supergravity (SUGRA) [U [2], super- weakly interacting massive fields are 
generally involved, such as the gravitino and moduli fields. They often cause cosmological 
problems, such as the moduli or Polonyi problem Hj, and gravitino problems [51 [5]. 
If they are generated in an early stage of the universe, they easily spoil the big-bang 
nucleosynthesis (BBN) when they decay, or tend to overdose the universe if stable. Hence, 
their abundances are severely restricted. Among the SUGRA particles, gravitinos are 
thermally produced in hot plasma. Moreover, it was pointed out that a large amount of 
gravitinos are generated via decays of coherent oscillations of heavy scalars, such as moduli 
and inflatons [H [91 [10] . Therefore, the cosmology constrains the reheating temperature 
and models of the scalars. 

The partial decay rate of a scalar into a pair of gravitinos depends on structures of the 
scalar and supersymmetry (SUSY) breaking sectors. The production rate has been studied 
in detail by Ref. [11| when the SUSY breaking is caused by a vacuum expectation value 
(VEV) of the F-term of the SUSY breaking field. In this case, the goldstino, which is the 
longitudinal component of the gravitino, is composed of the spinor component of the SUSY 
breaking chiral supermultiplet. In the literature which studied the gravitino production, 
D-term contributions to the scalar potential have been discarded. If the SUSY breaking is 
dominated by a D-term VEV, the goldstino mainly consists of a gaugino, and the above 
result can significantly change. In this paper, the direct gravitino production by the scalar 
decay will be studied when the D-term VEV is sizable. 

This paper is organized as follows. In Sec. the production rate of a pair of gravitinos 
will be provided, taking the D-term potential into account. The scalar mixing with the 
fields in the SUSY breaking sector will be investigated in the SUGRA framework. In 
Sec. El the result will be applied to a D-term SUSY breaking model, and its cosmological 
implication will be discussed. Sec. 0]is devoted to a summary and discussions. 

2 Direct Gravitino Production Rate 

In this section, we evaluate the partial decay rate of a scalar field such as the modulus and 
inflaton field into a pair of gravitinos. It is derived at the leading order of Planck-suppressed 
interactions in the framework of the SUGRA. Both the D-term and F-term SUSY breaking 
contributions are taken into account, and sources of the D-term contributions are clarified. 

2.1 Gravitino Production Rate in Mass-Eigenstate Basis 

The SUSY breaking sector is supposed to have nonzero VEVs of the D-term potentials for 
gauge symmetries which are not included in the Standard Model (SM) gauge groups. The 
gauge symmetries can be U(l) and/or non-Abelian symmetries. If the symmetry is U(l), 
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we assume that there is no genuine Fayet-Iliopoulos (FI) term, since it is difficult to embed 
the Fl-term in the SUGRA |12j . The VEVs of the D-term potentials are considered to be 
dynamically generated by the fields in the SUSY breaking sector, z l , which are charged 
under the gauge symmetries. See Sec. [3] for an example of such a model. 

In the gravitino production, the relevant fields are the scalar field, 4>, the SUSY breaking 
fields, z l , and the gravitino. It is assumed that the scalar <fi is much heavier than the 
gravitino due to the SUSY invariant mass0 The scalar is also assumed to be a singlet 
under the extra gauge symmetry, and couple with z % only through the terms suppressed by 
powers of the (reduced) Planck scale, Mp = 2.4 x 10 18 GeV, which generally exist in the 
SUGRA LagrangianJi 

The superpotential, Kahler potential and gauge kinetic function are represented ajf| 




where K^-L and W^^ are the interaction terms between <f> and z % whose Planck suppressions 
are displayed explicitly. The first term in the right-handed side of Uab corresponds to 
the gauge kinetic and 9 terms, and g is the gauge coupling constant of the extra gauge 
symmetry. Let us call the basis in Eq. ([1]) as the "model basis" in this paper. Note that 
the VEVs of (p and z % can induce kinetic mixings through higher dimensional terms in the 
Kahler potential. In the following, we assume that those VEVs are much smaller than 
the Planck scale. Otherwise, the following discussion applies after the fields are shifted to 
absorb the large VEVs. 

The decay rate is evaluated in the mass-eigenstate basis. The fields in the mass- 
eigenstate basis, X a , are related to those in the model basis, 

X a = A xa xa 5x a , (2) 

where 5x a = x a — {x a ) is the fluctuation around the VEV, and A is a matrix to diagonalize 
the mass matrix and to canonicalize the kinetic terms at the potential minimum. The 
model-basis fields, x a , consist of x a = <p a and its Hermitian conjugate tp a , with ip a = (p 

lr rhis is the case for the inflaton in many inflation models. In the case of modulus, if the modulus mass 
is < O(lOTeV), its decay occurs after the BBN starts and spoils the success of the standard cosmology. 
Thus, we assume ^> lOTeV > m 3 / 2 . 

2 Otherwise, (f> decays into the SUSY breaking sector more frequently, which worsens the cosmological 
gravitino problem. 

3 In this paper, we follow Ref. [2] for the convention. Derivatives with respect to fields are denoted by 
subscripts, e.g., G^, = dG/d<j>. Also, we omit symbols of VEV, (•••), if not otherwise specified. 
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and z l . On the other hand, the mass eigenstates are generally represented by real fields, 
X a = <!>/, Z R , and Zj, which primarily consist of the real and imaginary parts of 4> and 
z i , respectively. Hence, (A xa ^)* = A is satisfied. Since the interactions between cj> 
and z % are given by higher dimensional operators, A® R ' ! z i, A® R ' ! 2 i, A ZR ' I ( j ) , and A Zr ^^ are 
suppressed by the Planck scale. The matrix A will be evaluated later. 

The Lagrangian terms which are relevant for the tree-level decay of the scalar into a 
pair of gravitinos are found to be 
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(3) 



where ip a stands for (j) and z % in the model basis, while ip^ is gravitino, 7713/2 the gravitino 
mass, and G = K + ln|M^| 2 the total Kahler potential. Here and hereafter, the Planck 
unit, Mp = 1, is used if not otherwise specified. In terms of the mass-eigenstate basis, the 
above Lagrangian becomes 



where the coefficients are 

2i<Z(G tpa {A- x y° 



a (i) 



a (2) 



2K(G^( J 4- 1 )< $ 



(4) 



(5) 



where A -1 is the inverse matrix of A, and • • ) and £y(- • • ) represent the real and imag- 
inary parts, respectively. The omitted terms in Eq. Q include Z RI , which are irrelevant 
for the gravitino production by the decay. Note that the gravitino production rate 
will be derived at the tree level in this section. The rate can receive radiative corrections, 
which will be mentioned later. 

From the above interactions, the partial decay rate of the scalar into a pair of gravitinos 
is evaluated as [7] 



^(^R,i -> ^3/2^3/2) 



IG (eff) \ 2 m 5 
2887rm2 /2 



(6) 



where we have used m$ RI m 3 / 2 - The effective coupling constants are defined as 
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a (l) 
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y *R.. 



lfi a =(f>,Z i 



(7) 



The rate is apparently the same as that in the F-term SUSY breaking models [7J. The 
effective coupling constants (J7J) are governed by G^a , which is related to the F-term as 



F i = -e G/2 g fj G- i . 



(8) 
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where g l i is the inverse of the Kahler metric, gq = (Gjj). However, this does not mean 
that the D-term is unimportant. As we will see, magnitudes of F-terms are controlled 
by vacuum conditions, which are affected by the D-term potential. Moreover, the D-term 
contributes to the scalar mass matrix, which determines the mixings of scalar fields with 
the SUSY breaking fields. These D-term contributions to the effective coupling constants 
are studied in the following. 



2.2 Effective Coupling Constants 

In this subsection, the effective coupling constants of the gravitino pair production, , 

^ R,I 

are evaluated in the model basis within the SUGRA framework. Since we are interested 
in cosmological applications of the gravitino production, it is sufficient to evaluate them at 
the leading order with respect to the Planck suppressed couplings. Higher order corrections 
are safely neglected. The interactions are classified by powers of the inverse of the Planck 
scale. At the zeroth order, i.e., in the global SUSY limit, <p is secluded from the SUSY 
breaking sector by the assumptions. On the other hand, SUGRA corrections which include 
Planck suppressed interactions belong to higher orders of the perturbation. Turning on the 
corrections, the two sectors communicate with each other, and the gravitino production 
occurs. 

In Eq. Q, it is sufficient to evaluate G z i and (A~ 1 )^$, RI at the zeroth order, because 
and (A~ l ) z \ R r starts from the first order of the perturbation. Let us start from 
(A _1 )^$ HI . In the global SUSY limit, the mass eigenstates in the <j> sector are simply 
given by $r = A / g^/2(S(p + 



and 



/2(5(f) — 5(f>), with their masses 



1 



1 'v®-v& 



(9) 



where denotes the scalar potential in the global SUSY limit, and 

d 2 v 



dx a dx> 3 



(10) 



is the second derivative of the potential in the model basis. Here, we assumed vj^ = V~ 
for simplicity. It is straightforward to include the phase. Thus, the mixings of the scalar 
4> itself become (A~ 1 )^^ >R = 1 / ^J2g^ and (A -1 ) 9 ^ = i/y/2g^ at the zeroth order. The 
effective coupling constants ([7]) are approximated to be 



(eff) 
R 



(eff) 



VW 4 



-Ga) + G z i(A 



-l\z* 



(11) 
(12) 
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Next, the F-terms of the SUSY breaking fields, G z % ~ W z i/W, are evaluated at the zeroth 
order, namely by using the field VEVs in the global SUSY limit, as discussed above. This 
depends on models, and we will demonstrate it in a D-term SUSY breaking model in 
the next section. Although the SUSY breaking sector may be involved, this procedure is 
straightforward. 

The remaining task is to evaluate Ga and (A~ 1 ) Z \ R 7 at the first order. They include 
contributions from the D-term potential, which is represented by the Killing potential, Da- 
By the gauge invariance of the action, the Killing potential satisfies 

D A = iX A Gi = iX\K u (13) 

where X A = —ig lJ DAj is the holomorphic Killing vector, by which the gauge transforma- 
tion of a chiral superfield, <J? 4 , is defined as 5<fr l = A A X l A (&) with A A being the gauge- 
transformation parameters. The second equality in Eq. (|13j) is derived by the gauge invari- 
ance of the superpotential, 5W = Wi5<& % = 0. 

The F-term of <p is evaluated at the minimum of the scalar potential in the SUGRA. 
The conditions of the vanishing cosmological constant and the potential minimization are 

V = ^g 2 D A D A + e G (G i G l -3) = 0, (14) 

Vi = g 2 (^-h% i D A D B + D A D A )j + e G {G.G'G, - 2G t + G^ t G 3 ) = 0, (15) 

where V is the scalar potential in the SUGRA, h AB (h^) is (the inverse of) the real part 
of the gauge kinetic function, and ViGj = Gy — IyGfc is the covariant derivative with 

ij y ^ijl' 



Tfj = g^G^j. From Eq. (fPo|). Ga is obtained as 
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Y+G*, (16) 



"3/2 

where the gravitino mass is from m. 3 / 2 = (e^/ 2 ). In the expression, we have used [V^G^I = 

W^/W H ~ m0/m 3 / 2 > 1 and Gj < 0(1) from Eq. (fT4"|) . It is found that G<$, vanishes 

in the global SUSY limit, when <j> is secluded from the SUSY breaking sector in the limit 
and is a singlet under the extra gauge symmetry. 

Next, let us evaluate the mixing matrix (A~ 1 ) Z \ R r There are two sources of the 
mixing between cj> and z l . The first one is from the kinetic mixing, g z i^, which is induced 
by higher dimensional operators with the field VEVs. The kinetic terms are canonicalized 
by redefining the fields as 

9z i <j> ji iri-l\i 



+ ^^5z\ z' 1 = (C-yrfzi, (17) 



at the leading order of g z i^, where C is a matrix that canonicalizes the kinetic term of the 
SUSY breaking sector, 

CV'fei**C*, = S a . (18) 
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Here, it is sufficient to evaluate g zi ^ and by using the field VEVs in the global SUSY, 
when the first order perturbation is considered. The second source of the mixing comes 
from the mass term. By canonicalizing the kinetic terms, the mixings in the mass matrix 
become 



V, 



i 



:(CVtW, 



G 3 i V(f>zj ! 



(19) 



(20) 



at the leading order of g z i^. Explicit forms of V^ z i and V^ s j will be given later. Since 
the mixing terms are small, the mass matrix is diagonalized by means of the perturbation 
theory. At the leading order of the perturbation, (A~ 1 ) Z \ RI is obtained by combining the 
canonicalization and the diagonalization as 



i 



(A 



-l\z 1 



V^4 
i 



! 2 2\ — 1 

\™<s> R 9z ~ m z ) 



V&4> + Vsi<t> 



(mljgz 



m. 



) ] (VzH 



Vzi&> ! ■ 



(21) 
(22) 



where the indices <j> and z % are not in the primed basis of Eq. (|17p . but in the model one. 
Here, (m%g z 



m; 



is an inverse matrix, and the index z 3 runs over both z 3 and z 3 . The 



matrices, m 2 and g z , are defined in the model basis as 



nr. 



V (g) . v {g) . 

Z % Z3 z t z^ 

v {e) . U (g) . 

\ Z 1 Z3 Z % Z3 , 



9z 



9z i zi 

v g M i z j 



(23) 



Since the mixing V z j^ ± V^z belongs to at least the first order of the perturbation, the 
inverse matrix is evaluated in the global SUSY. 

The mixing terms V$ z i and V^ z i are obtained from the scalar potential of the SUGRA. 
By using Eqs. (|14[) and (|15p . the mass matrices are derived as 



Vr, 



e G (ViG k V 3 G k - R l]kI G k G T + 9i 



(24) 



+ g 2 Q (GiGj - 9i] ) D A D A - GiD A D A - - G ] D A D Ai + I (h%fi 3 + h%-G^ D A D B ^j 



<f ( h AB D Ai D Bl + hjFiD A D B3 + h AB 3 D A D Bl + D A D At] + h^h^h^D ' D 



■"ABj y 
R h CD h R 



A nB 



Vij=e G (2V i G i + G K V i V i G fc 



+ gf ( g ( G i G J ~ D D a - G{D D A j - GjD D A i 



(25) 



+^{h% i G j + hl Bj G i ) D A D B 



pfc 



-D A D Ak + —h ABk DD 



1 



+ g 2 ihjt B D A D Bt3 + h B B D A {D B j + hi B t D A D Bj + hj i B J D A D Bl + -h^ijD A D B ) , 
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with 



ViV,-G. 



j<^k 



Gijk 



(26) 



where the indices in a parenthesis are totally symmetrized. In each of Eqs. (|24p and (|25p . 
the first parenthesis is induced by the F-term VEVs, while the others are finite when the 
D-term contributes. At the first order, the mixing terms of the mass matrices which are 
(potentially) relevant for the cosmology are obtained as 



gWV+GfVgiGt + g zJ2 V^G z3 V- zk G- z , 
+9^^ '<j,G^\/ 'giGgj - R^izi^g* 3 * g z z G z iG z 
+ 9 2 { h R B DAcpD Bz -i + h R B ,pD A D m i + D D A(j> zi) , 
~ e ° (^'2 3 '^# + 2 9 z z Ggi^ z iG z ^ z kj g z z G z m 

+ 9 2 \ -g z3z G^ z kD A D Az j + hft B D A D B(t)Z i + hft B D A( /,D Bz i + h R lB ( f > D A D Bz 



(27) 



(28) 



Here, we have used \Ga\ <C 1, and kept only the terms which can be enhanced by the 
following reasons: (i) the modulus/inflaton is much heavier than the gravitino due to its 
SUSY mass, (V^G^I 3> 1, (ii) in the D-term SUSY breaking models, the fields in the 
SUSY breaking sector, z l , can have a large SUSY invariant mass, |V z iG zJ | 3> 1, and (iii) 
derivatives of D A with respect to z % can be enhanced, since the VEVs of z % are much smaller 
than Mp. 

In summary, the effective coupling constants are represented in the model basis as 

2 



(eff) 

R 



1 



(eff) 
I 



Ga, + G z i 

Ga + G ~i 



V&4> + Vzi$ 



{m^gz-miy 



(29) 
(30) 



where G^ is found in Eq. (|16fl . m 2 and g z are in Eq. (j23[) . and the mixing terms are 
in Eqs. (|19p and (|20p with ()27p and ([28]). They are independent of the matrix C in 
Eq. (fT7|) . The direct gravitino production rate is evaluated by substituting the above 
coupling constants into Eq. ([6]). 

The D-term contributions are found in the terms proportional to g 2 in Eqs. ()16p . (|27p . 
and (|28p . They are classified into two groups: i) contributions of </> to the gauge kinetic 
term such as h AB ~ (ft/Mp, and ii) those to the D-term potential, D A( j,, due to an effective 
charge which is induced by mixings with the SUSY breaking fields from higher dimensional 
operators. See Eq. (|13p for the expression of the D-term. Thus, all the D-term contributions 
vanish if the following two conditions are satisfied: 
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(I) 4> does not appear in the gauge kinetic function of the extra gauge symmetry, Kab$ = 
0, 

(II) (j) nas no kinetic mixings with the SUSY breaking fields which are charged under the 
extra gauge symmetry, = 0. 

If these conditions are satisfied, the gravitino production is only from the F-term contri- 
butions. (See App.[A]for a revision of the F-term contribution.) It is found that the rate 
decreases as the D-term VEV increases. This is because the VEV of the F-term potential 
satisfies the vanishing cosmological constant condition (|14p . If the D-term potential dom- 
inates the SUSY breaking, the F-term VEVs become suppressed. Then, the cosmological 
problem of the direct gravitino production can be relaxed. 

Finally, let us comment on radiative corrections to the pair gravitino production. The 
scalar field couples to a pair of gauginos via matter loops analogously to the anomaly, 
where the scalar interactions with the matter come from gravitational effects or higher 
dimensional operators [13} HOj . Since the gaugino in the SUSY breaking sector is a main 
component of the goldstino, the anomaly-induced decay is considered to contribute to the 
pair gravitino production. This decay induced by the gravitational anomalies works when 
the matter fields in the loop are lighter than the scalar. 

When the scalar 4> is heavier than the fields in the SUSY breaking sector, the latter 
fields other than the goldstino are also produced by the scalar decay through SUGRA 
interactions [21 [10]. The decay products in the SUSY breaking sector are considered 
to decay into lighter fields including the gravitino. This increases the abundance of the 
gravitino, and thus, worsens the cosmology. 



3 Example 

In this section, we apply the formulae in the previous section to a model of the D-term 
SUSY breaking. A fraction of the D-term SUSY breaking is defined as 

where the F-term SUSY breaking is represented by the F-term potential, Vf = e G G l Gi, 
and the D-term SUSY breaking is by the D-term potential, Vf> = (g 2 /2)D a Da- In the 
global SUSY, the F-term potential leads to Vf = g l:, WjWi. The denominator of the right- 
handed side of Eq. (|3ip is related to the gravitino mass through the cosmological constant 
condition (fH|) . i.e., 

V F + V D = g fj W 3 Wi + \g 2 D A D A = 3m 2 3/2 . (32) 

The fraction 5 takes a value of < 6 < 1, where a pure D-term SUSY breaking, 5 = 1, is 
excluded because of the gauge invariance condition, Da = D l A Gi [2]. 



S 



In Sec. 13. H a model of the SUSY breaking is introduced. The model is extended to 
include a scalar field (such as modulus or inflaton) , and the direct gravitino production rate 
is studied in Sec. 13.21 In particular, the conditions (I) and (II) obtained in the previous 
section are assumed to be satisfied. It will be explicitly shown that the rate is reduced 
compared to the F-term SUSY breaking models, as the D-term SUSY breaking effect 
increases. The cosmological implications are discussed in Sec. 13.31 

3.1 D-term SUSY Breaking Model 

Let us consider a model of the D-term SUSY breaking explored in Ref. [15J. The model 
has a U(l) gauge superfield and six chiral superfields, z-i,Zi,Z-xw,Zi/n,Zo and z' , with 
their U(l) charges denoted by the subscripts. The Kahler potential and the gauge kinetic 
function are minimal, while the superpotential is 

W{z % ) = Aiz (m 1 ~ N ZiZ J ? 1 / N - m 2 ) + \ 2 mz 1 z_ 1 + \?,ZqZi/ N z_i/ n , (33) 

where m determines the SUSY breaking scale. For A3 S> Ai S> A2 3> g, the field VEVs in 
the global SUSY limit are given by [15] 

( z _ 1 ) = (z 1/N ) = (z ) = (z' )=0, (34) 

\(z- 1/N )\^m(N*\t)^ , (35) 

(z 1 )^m N+1 (z_ 1/N y N , (36) 
where A2 = A2/5. At the vacuum, the scalar potential in the global SUSY is |15| 

V^^(N^(± + ^L). (37, 

Here, 1/./V 3 and 1/2N 2 in the parenthesis correspond to the F-term and D-term potentials, 
respectively. Thus, the model involves the D-term SUSY breaking effect as well as that 
of the F-term, and the fraction of the D-term SUSY breaking is approximately given by 
S ~ N/(N + 2). The D-term contributions increase as N becomes larger. 

In the following numerical calculations, we use more precise values of the field VEVs, 
which are obtained following the analysis in App. [Bl and the mass spectrum, which is 
obtained using the mass matrices listed in App. O The nonzero VEVs and mass spectrum 
are shown in Table Q] for some model points, as well as the SUSY breaking fraction 5. The 
numerical results in Table [1] agree with the approximate ones within ~ 20%. 

As shown in the table, all the fields in the SUSY breaking sector have masses scaled by 
m, which is related to the gravitino mass m 3/ / 2 by Eq. (|72[) . or approximately 

m ~ 6^N^g~^m 3 / 2 M P , (38) 

for large N, where we have used Eqs. ([32]) . ([35]) . and ([36]) . Note that there is no light scalar 
field (the so-called Polonyi field) in the SUSY breaking sector with a mass of order the 
gravitino mass. 
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A3, Ai, 












mass spectrum 




A2,9 


AT 


5 


<*> 




<*-l/JV> 


JSUSY sector 




PI 


1,10-!, 


1 


0.31 


0.46 




2.2 


3.7 x 10- 3 - 2.2 


3.2 x IO- 3 m 2 /M P 




io- 2 ,io~ 3 


10 


0.83 


8.3 x 10" 


-3 


1.6 


5.6 x 10~ 4 - 12 


1.2 x 10- 4 m 2 /M P 






100 


0.97 


1.4 x 10" 


-4 


1.1 


1.3 x 10" 4 - 7.0 x 10 2 


5.0 x 10" 6 m 2 /M P 


P2 


1 4" 1 


1 


0.27 


0.61 




1.6 


4.0 x 10~ 2 - 1.6 


2.6 x 10" 2 m 2 /Mp 




4- 2 ,4- 3 


10 


0.83 


1.8 x 10" 


-2 


1.5 


8.1 x 10- 3 - 14 


1.6 x lQ- 3 m 2 /M P 






100 


0.98 


2.9 x 10" 


-4 


1.1 


1.7 x 1Q- 3 - 8.5 x 10 2 


7.6 x 10 _5 m 2 /M P 



Table 1: The fraction of the D-term SUSY breaking 8, the field VEVs, and the mass 
spectrum for N = 1, 10, and 100. Here, two sets of the parameters, PI and P2, are chosen. 
See Apps. [B] and [C] for details of the analysis. The field VEVs and the masses are written 
in units of m except for the gravitino mass m^^- 



3.2 Gravitino Production Rate 

The model in the previous subsection is extended to include the scalar <f>. For simplicity, 
it is assumed that the Kahler potential and gauge kinetic function are minimal, 

K = \(j)\ 2 + \z \ 2 + |zq| 2 + \zi\ 2 + |z_i| 2 + I^i/tv) 2 + |^-i/Af| 2 , h u{ i) = l. (39) 

It is also assumed that the mixings between eft and z % are absent in the superpotential, 
wj^l((/), z l ) = 0, and the total superpotential is given by 

W = W(z i ) + W(<f>), (40) 

with W(z l ) given by Eq. ([33]) . It is found that the conditions (I) and (II) discussed in the 
previous section are satisfied. 

Let us now calculate the partial decay rate of the scalar <f> into gravitinos, or the effective 
coupling constants in Eqs. ([29]) and ([30]) . following the discussion in the previous section. 
They are determined by and the mixings between the scalar and the SUSY breaking 
fields. In the present model, from Eq. ([To]) . becomes 

G \.i ^7 J.G ~i 

G *-^f' (41) 



and the mixings in Eqs. (02]), flUJ), ([27]), and ([28]) are 

V+gi = V+gi ~ e G (V^V^ + Va« 2 iV 2 i^V) , (42) 

%* = V - °- ( 43 ) 
In Eqs. (jHJ) and ([4"2]) . the factors e G , V^G^, and V z iG z j are zeroth order in the perturba- 
tion, which are given by 

V, e ,-«5*-^. (44) 
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On the other hand, the factor V ' z iG^ = V ^G z % is suppressed by the Planck scale. Due to 
the absence of the <f> — z l mixings in K and W in the present model, it is given by 



(K^ - G 4 



w 



K, 



' w 



(45) 



w w 

Here, in the last equality, we have used G^ ~ x 0(m 3 / 2 / m 4>) ^ K<f>, which can be 
shown by using Eq. ([IT]) . V^G^ ~ m^/m^^i an d G z i ~ W z i/W < 1. See also App. [A] 
The remaining task is to evaluate the F-term VEVs of the SUSY breaking fields, W z %, in 
the global SUSY limit. From Eqs. ([33]) - ([36]) . the finite F-terms are 



X 1 (m 1 ~ iV ( Xl ) (z_ 1/N ) 



N 



m 



ir 



2-1 



A 2 m (zi) 



(46) 



it is found that W z _ x dominates the F-term potential]^ The other 



From Eqs. 

F-terms are zero in the global SUSY. From Eqs. ([13]) and ([15]) . it is found that V^G^ is 
dominated by V 'd>G z _ 1 as 



m - 



3/2 

where (W) — rn S fy is used. Thus, the couplings are dominated by z % 
the above equations altogether, one obtains 



,(eff) 



1 



m 



I)" 



Z_12_1 



^3/2 



G z _^7 JG Z 



(47) 
Putting 

(48) 



For the purpose of the cosmological application in the next subsection, we concentrate on 
the case in which the scalar (j) is lighter than the fields of the SUSY breaking sector z % . 
Then, the effective coupling constants are reduced to 





2 







in 



3/2 



G Z _ 1 V 4> G Z _ 1 



Xlm 2 



36 ™ 3 /2 ,.,2 



rn 



(49) 



In the last equality we have used Eqs. ()35p and ([36]) . and N 3> 1, where the effective 
couplings become insensitive to the model parameters of the SUSY breaking sector except 
for N. 

In Fig. [JJ we show the effective coupling constants as a function of the fraction of the 
D-term SUSY breaking 6, which is varied by changing N of the model. Here, we have used 
the full formulae ([29]) and (|30]> with the numerically obtained field VEVs|j The mass of 
the scalar is set to be smaller than the masses of the fields in the SUSY breaking sector. 
For the sake of comparison, we have also shown, by the star sign in the figure, the case 



4 Although W ZQ seems to vanish when the VEVs in Eqs. I|34[) _ (|36[) are naively applied, it is finite as 
long as the D-term is nonzero, because the right-handed side of W ZQ is analytically related to the D-term 
VEVs via the stationary condition of the scalar potential (|71[1 . The gaugino masses of the SUSY SM can 
be given by W ZQ , or by the anomaly mediation. 

5 We have also used Eqs. (|44[1 and (|45(l to evaluate the derivatives of W(4>)- 
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(eff) 
Ft 



(eff) 



on the 
10 12 GeV 



Figure 1: Dependence of the effective coupling squared 

fraction of D-term SUSY breaking 5. = 10 3 GeV, = 10 6 GeV, and (cj)} 

are used. The lines for the parameter sets PI and P2 in Table Q] coincide. The fraction 
5 has been extrapolated to N < 1, which is drawn by the dashed line. The star mark 
represents the case of the pure F-term SUSY breaking. 



where the SUSY is broken only by a single F-term VEV, assuming the minimal Kahler 
potential: 



,(cff) 



G Z X7 aGz 



3m 



3/2 



<</>> 



(pure F-term), 



(50) 



where we have used | G z | ~ y/3. It is found that the effective coupling constants diminishes 
as the D-term SUSY breaking contribution increases. 



3.3 Cosmological Implications 

Finally, we estimate the total gravitino abundance. The coherent oscillation of the scalar 
4> is assumed to dominate the energy density of the universe and then decay to create the 
thermal bath. There are two types of gravitino productions: one is the thermal production, 
where the gravitinos are produced through scattering processes in the thermal bath, and 
the other is the direct production, which is the main subject of this paper. In the fol- 
lowing, we focus on the parameter region where the scalar is lighter than the fields in the 
SUSY breaking sector except for the gravitino. Otherwise, the anomaly-induced gravitino 
production and the direct production of the SUSY breaking fields which are mentioned in 
the last section can spoil the success of the standard cosmology. 

First of all, let us consider the thermal production of the gravitino. The gravitino yield 
is defined as Y s / 2 = n 3 /2/s, where n 3 / 2 is the number density of the gravitino, and s is the 
entropy density. The thermal gravitino yield is approximately given by [6] 



-v^thermal 
r 3/2 — 



1.3 x 1(T 14 + 



x 10" 




1 10 8 GeV 



(51) 
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where Tr = (7t 2 (7*(Tr)/90) 1//4 y/T^ is the reheating temperature after the (f> decay, with 
the total decay rate T^, and m\/2 is the unified gaugino mass at the grand unified theory 
(GUT) scalefl 

Next, let us consider the direct production of the gravitinos. The direct production of 
a pair of gravitinos provides 



Y 



decay 



3/2 



477U 



x 2B 



3/2 



90 



192vr \^g,(T R ) 



m. 



in 



3/2 



Tr 



(52) 



Here, m$ R ~ m<$, 1 = and 



(eff) 

R 



(eff) 
I 



are used. Substituting Eq. (jl9j) 



into Eq. (|52j) . the yield is reduced to be 

90 



y: 



decay 



1 



3/2 



192tt \ir 2 g*(T R 



~ 8.2 x 10 



-15 




(53) 



for large N, where ^(Tr) = 228.75 is used. It is found that the yield is proportional to 
m 2 , ((j)) 2 , and Tf, 1 . In particular, it becomes independent of 7743/2 because (zi) oc m oc 
-^/ m 3/2- It is stressed that the abundance decreases as N increases, namely the D-term 
contributions dominate (See also Fig. [T|). 

The cosmological constraint on the gravitino abundance depends on the gravitino mass. 
When the gravitino is unstable, the most stringent bound comes from the BBN. The precise 
bound depends on details of the mass spectrum including the SUSY SM. Here, we adopt 
the constraint for the model point in Case 2 in Ref. [6J0 where the bounds are given 
by Y 3 / 2 < 2 x 10~ 16 and Y 3 / 2 < 2 x 10~ 12 for the gravitino mass of m 3 / 2 = ITeV and 
30TeV, respectively. Thus, from Eqs. ([5T]) and (|53|) . the allowed region of the reheating 
temperature is 



1 



3 x 10 D GeV x 



(- 



m„ 



N 2 \ 10 12 GeV / Vl0 12 GeV/ 



Y < T R < 2 x 10 6 GeV (m 3/2 



5 x 10 2 GeV x 



1 



N 2 V 10 12 GeV 



(0) 



in,, 



10 12 GeV/ 



) 2 < T R < 1 x 10 10 GeV (m 3/2 



= ITeV), 
(54) 

30TeV), 
(55) 



where the upper bound comes from the thermal gravitino abundance, while the lower bound 
is from the direct production. This is because the thermal abundance is proportional to 
Tr, while the abundance from the latter is to 1/Tr. It is noticed that the window tends to 



6 The definition of Tr is different from the Ref. [6\ Logarithmic corrections are omitted in Eq. (|51|) . 
7 This model point is excluded by the recent SUSY search at the LHC, but we adopt it just for illustration. 
The main conclusion does not depend much on the details of the mass spectrum. 
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10 4 10 5 10 6 10 7 10 8 10 9 
rrifi [GeV] 

Figure 2: Contours of the cosmological constraints on the - (0) plane from the thermal 
and direct productions of the gravitinos. The upper-right regions are excluded. Several 
values of the D-term SUSY breaking contributions are taken: 5 = (the pure F-term; see 
Fig. Q] and the text), 0.27 (N = 1), and 0.98 (N = 100). The dashed lines correspond to 
777,3/2 = lOGeV, and the solid ones correspond to 7773/2 = ITeV. The former lines end in 
the figure because is restricted to be smaller than the smallest mass of the fields in 
the SUSY breaking sector. The reheating temperature is provided by = m^/8-K. The 
thermal abundance is subdominant in this figure. 

close as the energy scale of the scalar model is higher. In order to avoid the cosmological 
bound, the D-term contributions must be enhanced by increasing N. For instance, the 
window opens only when N S> 1 for 7773/2 = ITeV, ((/)) = 10 16 GeV, and = 10 9 GeV. 

When the gravitino is the lightest among the SUSY particles, the gravitino becomes 
stable. Since the massive gravitino contributes to the dark matter abundance, the gravitino 
abundance cannot exceed the cold dark matter abundance as 

^3/2^3/2 < ^dm < 4.4 x l(T 10 GeV, (56) 

where p c is the critical density, and Odm^ 2 < 0.12 at 2a [T6| is used. For the light stable 
gravitino, the bound becomes 

(57) 

where the gaugino mass is taken to be 600GeV at the GUT scale, which affects the upper 
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10 4 10 5 10 6 10 7 10 8 10 9 
rrifi [GeV] 

Figure 3: Same as Fig. [2] except that the reheating temperature is optimized such that 
the total gravitino abundance is minimized. The bound on m<f, and {(/)) is considered to be 
most conservative. The gravitino mass is set to be m 3/ / 2 = ITeV. There are no bounds in 
the case of 7713/2 = lOGeV in this parameter space. 

bound via the thermal gravitino abundance. Again, the lower bound becomes milder as 
the D-term contributions dominate. 

Finally, let us show how the D-term SUSY breaking contributions enlarge the allowed 
parameter space spanned by and (<f>). In Figs. [2]and[31 contours of the cosmological 
bounds are drawn for various S in the - {(ft) plane. The region of lower and {4>) is 
allowed. In the figures, is restricted to be smaller than the smallest mass of the fields in 
the SUSY breaking sector, for simplicity. It can be seen that as the D-term SUSY breaking 
contributions increase, i.e., as 5 becomes close to unity, the allowed region becomes wider. 

The bound depends on the reheating temperature as seen in Eqs. (|5ip and (|53p . In 
Fig. [21 the temperature is estimated by assuming that the total decay rate is governed by a 
Planck-suppressed operator, T((j) — > all) = m^/87r. The thermal production is subdominant 
in the figure. On the other hand, the reheating temperature is taken to be most conservative 
in Fig. [3l Since the thermal gravitino abundance is proportional to Tr, and that from the 
direct production is to its inverse, the total gravitino abundance becomes the minimum 
value by tuning Tr. Therefore, the bound in Fig. [3] is the most conservative one, and the 
model is refuted if it is excluded in Fig. [31 unless there is an entropy production. 
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4 Summary 



In this paper, we studied the scalar decay into gravitinos in the presence of the D-term 
SUSY breaking and multiple SUSY breaking fields, taking into account the mixing between 
the scalar and the fields in the SUSY breaking sector. We obtained the general formulae 
of the effective coupling constants for the gravitino pair production, (|29p and ()30p . It is 
found that the gravitino production is suppressed if the conditions (I) and (II) in Sec. 12.21 
are satisfied. 

As an example, we applied the formulae to the D-term SUSY breaking model in 
Ref. [15], and showed that the D-term SUSY breaking suppresses the gravitino production. 
Moreover, an interesting feature of the model is that all the fields in the SUSY breaking 
sector have masses of 0{- s J M p m^/2)- This is very different from F-term SUSY breaking 
models, in which the SUSY breaking field tends to have a mass of 0(m 3/ / 2 ); causing the 
Polonyi problem [3]. This may be a generic feature of D-term SUSY breaking models. (See 
also a recent work |17j.) These observations can open a new window for model building 
respecting cosmology. 
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A Gravitino Production by F-term SUSY Breaking 

In this appendix, the direct pair gravitino production by the scalar decay is revisited 
when the SUSY breaking is caused by the F-term VEV. The derivation of the decay rate 
explored in this paper is somehow different from the previous one [11]. In particular, the 
SUSY breaking sector is generalized, and the source of the gravitational corrections are 
clarified. In this section, the previous result will be reproduced. Moreover, the contribution 
which was overlooked will be pointed out. 

It is assumed here that the SUSY breaking sector is composed of a single chiral su- 
permultiplet, z, and there is no extra gauge symmetry for the sake of comparison with 
the result in Ref. [TT]. Also, the kinetic terms are supposed to be canonicalized in the 
global SUSY limit. In contrast to the D-term SUSY breaking, the SUSY invariant mass 
of z is at most comparable to the gravitino mass, V Z G 2 ~ 1, which was neglected in Ref. 
[11] . However, the scalar mass of the SUSY breaking field, m z , can be lifted by a higher 
dimensional operator, 5K = — |z| 4 /A 2 with A <C Mp. Under the assumptions, the effective 
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couplings (|29|) and ([30]) become 



,(eff) 



Gs + G 5 



m« — ml 

<2R,I 2 



(58) 



where and the mixing terms are simplified to be 

1 



Gr! 



v m 



VzG. 



-G z VjG S + QrbzG, 



(g) 



e G {V A G 6 V Z G^ - gtfV+G+VgGg - R rzz - Z G- Z G Z ) - g&V^ , 



9ziK 



(g) 



eG Gs zz GssGz - g z<b V?^ . 



(g) 



The mass of the scalar, 0, is 



66 4"P 



(g) 



m 



(59) 

(60) 
(61) 

(62) 
\. which is 



In particular, is dominated by the superpotential term, 

related to the SUGRA term as e G ^ 2 V ( j ) G ( f ) ~ W^. On the other hand, the last term 



is small, which originates, e.g., in the SUSY breaking effect, (= AmJ) ~ m 3 /2 m 6- 
The effective coupling is compared to that in Ref. [11] by considering the following, 



2 _1 

= 2 



,(eff) 



+ 



,(eff) 



mz 



m 6 ~ 771-2 



-a 



G Z S7 Z G^ 
ViGi 



zG z 



vnh — ml 



+ 



e {g z j?J 6G&V Z G Z + RdzzzG z G z ) 

2 



G 



2 2 9 

ml — m. 



G (hzsG \aG ' 2 



(63) 



where the corrections proportional to AmJ from (w-| fi 7 — m^) -1 partially cancel those from 
Vj^, and the remaining corrections are subleading and omitted. This result reproduces the 
gravitino production rate in the F-term SUSY breaking models |11| . 

Among the contributions in Eq. ([63]), V Z G& is represented in terms of K and W as 



V Z G, 



K .Wz ± _W z W 6 _ i 
^ W WW L *r' i > 



where W&/W can be replaced by 



w *-r k 



-K<k + g<bzG z , 



(64) 



(65) 
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up to a correction of order m 3 /2/ m 4>- This is because 
W<p m 3/2 



rriA 



G z V^Gz + g<j,zG z — K§ 

G z (k, $ + V ^-^-- 7 - Y'G, I +< lll .G, - ft',, (GG) 



W m. 

™3/2„ ( 11 IF- IV, 



and hence 



1 + ° ( ^t)) w = + ^ ! — - r, a, ) . (07) 




J J W 
Thus, Eq. (fM|) becomes 

IV^ w 7 w, 

V,G ~ tf^ + + - 94,-zG^ - T^Gz. (68) 

Barring accidental cancellations, V ' z Ga, ~ {4>) is obtained especially from the third term 
in the right-handed side. In Eq. (|63p . the terms in the first parenthesis of the most right- 
handed side, which include the V ' Z G^ term, are found to be suppressed if the scalar mass, 
m?, is much larger than that of the SUSY breaking field, m^. They originate in and the 
mixing V^, which cancel each other, as pointed out in Ref. |18| . However, this cancellation 
does not work for < m z , and the direct gravitino production takes place generally. On 
the other hand, the other terms in Eq. (|63p depend on the higher dimensional operators. 

The term proportional to g z zS/ '^G^V 1 z,Gz in Eq. ([55]) is missed in Ref. [TT| because 
V s Gz was supposed to be zero. If it is finite in the global SUSY limit, the contribution 
can be relevant when m<h is larger than m z - 



B Potential analysis 

In this appendix, the scalar potential of the model in Sec. [3] is analyzed. First of all, the 
parameters are taken to be real and positive by redefinition of the fields. In the global SUSY 
framework, the phases of fields are minimized, except for the would-be Nambu-Goldstone 
(NG) boson. Next, the VEVs of zq and z' can be represented in terms of those of the other 
fields, since the potential is quadratic with respect to zq and z' . Substituting them into 
the potential, it involves four real variables, z\, Z-\, z 1 / N , and Z-\m. It has quartic and 
quadratic terms with respect to z~\. The coefficient of the z^_ x term is positive, and the 
coefficient of the z 2 _ x term is also positive if the following condition is fulfilled, 

N^z^Xl 



■2 



9 Z (%i> 2 + ^ ({zi/n? - {z-i/n?) ) > 0. (69) 



If this is fulfilled, the VEV of z~\ vanishes. Similarly, the VEV of Z\m vanishes if the 
following condition is satisfied, 

A (z-i/n) 2 + ^ ((zi) 2 - i (z-i/n) 2 ) > 0. (70) 
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The VEVs of zq and z' also vanish if these conditions are satisfied. Then the potential 
becomes a function of z\ and z_i/jv> 



(71) 



Finally, we obtained the VEVs of z\ and z_i/at (in units of m) numerically by minimizing 
the potential, and confirmed the two consistency conditions (f69|) and ([TO]) are satisfied. 

Once the field VEVs are obtained, the SUSY breaking scale can be written in terms of 
the gravitino mass, by using Eq. (|14p or ([3"2j) : 

(2 2\ -1 / 4 

X\{h^ VN -l) 2 + \lzl + 9 -(^l-^z\ /N ^ J ^m 3/2 M P , (72) 

where and £_i/jv denote the absolute value of each VEV in units of m. 

C Mass-matrix components at the vacuum 



In this appendix, we explicitly write the mass matrix of the model used in Sec. [3] at the 
vacuum in the global SUSY. In the following, all the parameters and VEVs have been 
transformed to real and positive variables. In this appendix, we take m = 1, and use a 



shorthand notation for field derivatives like V. 



V 



Z % Z3 ' 



C.l Mass- matrix components of the scalars of the SUSY breaking sector 

The scalar mass matrix can be written in a block-diagonal form, 

/ A 



where 



B 
\ C 



A 



B 



C 



( U (g J 
1,1 

rig) 




V 



V- 



v l,l 
(g) 

1/N,1 

yjg) 
1,1 
(g) 



V, 



V 



(g) 

1,-l/N 

(g) 

l/N-l/N 



(g) 



1,-l/N 



V- 



(g) 



1/JV.l -1/N,-1/N ) 
V (S)_ y(g) 



1,-1 

u (g -L 

0,-1 



if 

0,0 



r(g) 



with c = U^ B; — and c' = VI ,. 

1/N,1/N 0',0' 







(73) 



(74) 



(75) 



(76) 
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The nonzero components are given below, 



V^=Xl + Xlz^ /N + 9 ^(2zl-^_ 1/N ). (77) 



V (S) _Ar2 A 2 2 2(JV-1) 9 2 ( 2 2 2 \ /yo^ 



V 



(g) -NXl Zl ^ - f^-i/iv, (79) 



l-l/JV 

= 5 2 , 2 , (80) 



y%,-y N = N ( N - ^W-i/n (*i*?i/N - l) + I^V' (81) 

=^?*5~/Jr (**V - l) - W, (82) 

^=Ai-^ 2 (^-^V), (83) 

= A i^V, (85) 
« 2 / 1 

v2„2 i 5 / 2 1 „2 



^vT7v=^l/* + ^ Ul-^-W , (86) 



/(g) _\2 7 2 



K% =Xlz\ IN . (87) 



The mass matrix has one zero eigenvalue, which corresponds to the would-be NG boson of 
the U(l) gauge symmetry. 

C.2 Mass-matrix components of the fermions 

The nonzero components of the mass matrix components Mjj are listed below. The sub- 
script "U(l)" denotes the gaugino. 

M_ M =A 2 , (88) 

Mi, =\iz N l/N , (89) 

M_ 1/JV , =NX 1 z 1 z N y 1 N , (90) 

M liU{1) = - V2ig Zl , (91) 
\/2 

M -i/A^(i) = jfigz-i/N, ( 92 ) 

M 1/JV)0 ' =A 3 ^i/tv- (93) 

The mass matrix has one zero eigenvalue, which corresponds to the would-be goldstino of 
the local SUSY. 
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